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Abstract 

The behavior of the two-particle Green's function in QED is analyzed in the limit 
when one of the particles becomes infinitely massive. It is found that the dependences of 
the Green's function on the relative times of the ingoing and outgoing particles factorize 
and that the bound state spectrum is the same as that of the Dirac equation with the 
static potential created by the heavy particle. The Bethe-Salpeter wave function is also 
determined in terms of the Dirac wave function. The present result excludes the existence, 
in the above limit, of abnormal solutions due to relative time excitations as predicted by 
the Bethe-Salpeter equation in the ladder approximation. 

PACS numbers: 11. 10. St, 12.20. Ds. 

Keywords: Relativistic bound states, Bethe-Salpeter equation. Quantum Electrody- 
namics. 



^E-mail: jallouli@ipno.in2p3.fr . 
^E-mail: sazdjian@ipno.in2p3.fr . 

^ Unite de Recherche des Universites Paris 11 et Paris 6 associee au CNRS. 



1 



1 Introduction 



It is generally admitted that the infinite mass limit of the two-particle Green's function in 
QED (or in similar abelian theories) yields for the lighter particle the Dirac equation in the 
presence of the static potential created by the heavy one. The corresponding derivations 
are essentially based on mass-shell |T|, ^ |^ or quasi-mass-shell ^, ^ conditions. On the 
other hand, it is known that the Bethe-Salpeter equation in the ladder approximation 
leads to the appearance of "abnormal" solutions in the spectrum, due to the relative time 
excitations [|^, ^, ^, |T0[. Taking the infinite mass limit (for one of the particles) does not 
remove from the spectrum the above solutions and therefore the corresponding limiting 
equation is not equivalent to the Dirac equation with a static potential. It is noteworthy 
to mention that in the static limit, when both particle masses go to infinity, the abnormal 
solutions still persist in the ladder approximation Jll] . 



In a recent work, we showed that in the static limit the two-particle Green's function 
can exactly be calculated and the corresponding bound state spectrum does not display 



any abnormal solution |]T2[. This result strongly suggests that the appearance of the 
abnormal solutions is rather a consequence of the ladder approximation (with covariant 
propagators) of the Bethe-Salpeter equation. It is then natural to expect that also the 
same conclusion might hold in the one-particle limit. 

The purpose of the present article is to study the (one-particle) infinite mass limit 
of the two-particle Green's function without making any mass-shell or quasi-mass- shell 
assumptions on the heavy particle and therefore not constraining the relative time evo- 
lution of the system. We shall show that the dependences of the Green's function on the 
relative times of the ingoing and outgoing particles factorize and that its exact spectrum 
in the above limit (in the absence of radiative corrections) is given by that of the Dirac 
equation of the lighter particle in the presence of the static potential. The relative time 
variables, without being absent from the expressions of the Green's function and of the 
Bethe-Salpeter wave function do not induce energy excitations, but rather ensure by their 
presence the correct spectral properties as dictated by quantum field theory ^. This 
result also also allows the determination of the Bethe-Salpeter wave function in terms of 
the Dirac wave function. 

While in the above calculations radiative corrections of the lighter particle are ne- 
glected (or partly taken into account in covariant gauges), this approximation does not 
seem crucial for the derivation of the factorization property of the Green's function in the 
relative time variables, since this feature is rather related to the heavy particle proper- 
ties. By appropriately generalizing the static potential so as to include the full radiative 
corrections of the lighter particle, one should still be able to prove the main qualitative 
results as described above. 



2 



The QED bound state spectrum being gauge invariant it is actually sufficient 

to prove the above properties in one particular gauge. We shall choose the Coulomb gauge 
for a detailed derivation, but shall also briefly sketch the line of calculation with linear 
covariant gauges. 



2 The two-particle Green's function with an infinite 
mass 

Our starting point is the QED lagrangian density for two fermions with masses mi and 
1712 and charges ei = 62 = — e: 

C = ipi{i'y.d — nil + (^l-^)^^! + V'2(n-"9 — m2 + ej.A)ip2 

where £}.{. is the gauge fixing part of the lagrangian density. The two-particle Green's 
function is defined by the functional integral: 



^ < M^i)Mx2)M^2)M<)e''^^'^'^^ >, (2) 

where S = J d^xC{^) is the action. We shall be interested by the bound state spectrum 
of the system fermion 1-antifermion 2. We consider particle 2 as the heavy one and 
take the limit of large m2. This limit can be studied by treating perturbatively the 
(covariant) kinetic energy term of particle 2. Such a situation has often been met in 
heavy quark effective field theories [|T3|. A way of implementing this limit is to make 
successive transformations on the fermion 2 field, the ffist of which being |T^ : 



i/j2 exp [{i-fid' + e-fiA')/{2m2)) ^2, ^2 ^ ^2 exp ((-^7* +67^ A0/(2m2)) . (3) 



This transformation has been shown to have unit jacobian UMf. It eliminates at order 
m2 the (covariant) kinetic energy part of the fermion 2 lagrangian density and the corre- 
sponding free propagator becomes the static one: 

G2o(4 - X2) = \{l + 7o)e-^'"^(^'"-^")^(x'2° - xl)5\^^ - x^) 

+i(l - 7o)e-^'"^(^^^'")^(x° - x'^)6\^2 - A)- (4) 

The antifermion contribution is selected when x\ > xf and we shall henceforth consider 
only this case (70 has then the eigenvalue —1). 
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Integration of the functional integral with respect to the photon field amounts to 
replacing the action S by an effective one S = J d'^xd'^yC{x,y), with 

C{x, y) = (^'1(27. 9 - mi)^i + ii2{ilod^ - m2)ip2)^ S^{x - y) 

+ ^ + J2,{x)) D^^{x - y) [jUv) + J2M) , (5) 

where j^/. = ^a7/.V'a, « = 1,2, is the current of particle a, with j2/i = fl';,°j2o, and D^^ is 
the photon propagator. We shall neglect in the following tha radiative corrections and 
therefore shall maintain in the interaction part of C the crossed terms ji-Dj2 + J2-Dji 
only. (Actually, the radiative corrections relative to the heavy particle 2 can explicitly be 
calculated |12| and the latter approximation concerns rather the lighter particle 1.) 

Use of Wick's theorem with respect to particle 2 together with formula and con- 
dition X2 > X2 yields for the Green's function (H) the expression: 

G{xi,X2,x[,X2) = G20{X2 - X2) < ^1 (xi (x'l) 

X expjz J d^yiij^{yi)(^ijiA - mi - ie^-f'^Af,{yi,X2,X2)^tpi{yi)^ >, 

(6) 

Af,{yi,X2,X2) = / dy^D^Qiyl - ?/°, yi - Xa). (7) 
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3 Coulomb gauge calculation of the energy spectrum 

In the Coulomb gauge, is: 



4tt |yi -X2I 



x',-y",)e{yi'4l- 



The Green's function (|^ can be evaluated by considering the differential equations it 
satisfies. These are: 



G{xi,X2,x[,x'2 



= G2g{x'2 - X2)i5^{xi - x'l), (2) 
- Z710 4o -^7i 4 -nil - lioV{r')9{xl - x'°)0(a;'° - x'^) 

= G2o{x2 - X2)i6^{xi - x[), (3) 



a 



r' = Ix; - x'. 



(4) 



where 

V{r) = , « = 1— , r=|Xi-X2|, 1 -1^1-^2 

r Att 

Before solving equations (^-(^ we shall introduce a particular representation for the 
one-particle Green's function which will be useful throughout the resolution procedure. 
Let us first consider the free Green's function of particle 1. Defining in momentum space 
the functions 

^o(Pi) = 7o(7iPl + mi) 

and the projectors 



ml + pI 



one can write the free Green's function in the form: 

i . i 



(5) 
(6) 



a 
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Pio - ho + ie 



Ao+70 + 



Pio - ho-ie 



Ao-7o- 



(7) 



It takes in x-space the form: 



Gioixi - x[) = e-^'^U e{x\ - ) + - Aho^<"'''", 



where now ho and Hq are the x-space expressions of the function /io(Pi): 

^0 = --ilol'dii + mi7o, h'^ = -i-foY d[i +mi7o. 



(9) 



By replacing, through the closure relation, the delta-function above by a complete set of 
spinor wave functions and identifying the arguments of the exponential functions with 
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the energy eigenvalues of these functions, one deduces that for x'J > xf they correspond 
to the positive energy solutions of the free Dirac equation and for < x'^ to the the 
negative energy ones, thus exhausting the spectrum of states saturating the free Green's 
function. 

The above construction can also be generalized to the case of the one-particle Green's 
function in the presence of a static vector potential V{r). The Green's function then 
satisfies the differential equation 

(^i'ji.di - mi - 7ioV(r)^Gi(xi,x'i,X2) = W^(xi - x'l) (10) 

and its conjugate one. Since V is static, the formal solution of these equations can be 
obtained from expression (|]) by a shift of the operators and /iq by the factor V. 
Defining: 

h = -t-foi'du + mi7o + V{r), h' = -i^oY d[i +mi7o + V{r'), (11) 

one obtains: 

Gi(xi,x;,x,) = e-^''^i(e(x? - x'°) + -^^'(^1 " xi)7oe'^>'. (12) 

[a//? is an appropriate generalization of \fh^.] We also notice the following property of 
the operators h and h': 

M=^(xi - x;)7o = (5=^(xi - x;)7oh', (13) 

which allows us to convert, when necessary for symmetry reasons, the operators h and h' 
into each other. Again the use of the closure relation through the delta-function allows us 
to identify the spectrum os states saturating the Green's function with the set of solutions 
of the Dirac equation in the presence of the static potential. 

We are now in a position to construct the solutions of Eqs. (0)-(^. We first rewrite 
them in the form: 

G = G2oi6\xi - x[hio, (14) 

= G2oiS^{xi - x[)-fio. (15) 

After making the change of function G = ^-i-A^Q^'^^i^' ^ one integrates for G, by imposing 
the boundary condition that in the limit x'2 — > — 00 and x^ 00, with x'^ and x\ fixed, G 



idw -h + V[l- e^x'i - x\)e{xi - x'^ 



G 



- i d[o -h' + V\il- e(a;° - x'^)e{x'l 



/o 



^2 , 



6 



tends, up to multiplicative factors, to the static potential case solution (JT^), as is evident 
from Eqs. (0)-(^. One finds the solution: 

G = G2o(4 - X2)e-'''>T exp i^i rfyje**'? V(r)(l - - ?/?)^(?/? - 4°))e-^^?'^| 

xl{ei4 - + ^y'(^^ - ^hioe^'"'"''' , (16) 

where T designates the chronological product. This expression can be further simplified 
for the situation that interests us. Sticking for definiteness to the case X2 and greater 
than x'2 and xf, one observes that the domain of integration [0:2°, 0:2] does not contribute 
to the above integral. One therefore obtains two disjoint domains of integrations, which 
allow the factorization of the T-product. Using then the property (|13|), the expression of 
G becomes: 

G = ^20(4 - a;2)e-^^?''Texp dy^^e'y>V{r){l - - - 4°))e-^^?'^' 



x^fe(x? - ) + ^ )()"(^xi - Xij7io 



2V ^ ^ ^' ^h^J 
xTexp \i J J rfyiV^?'^V(r')(l - - y'MVi ' 4°))e-^^'|. (17) 

Introducing the relative variables x = X1—X2 and 2 , and making in the integrals 

appropriate changes of variables, and also using the property exp(e'^y4e~'^) = e^e'^e'^ , 
one finally obtains for G the expression: 

G = G2o{x2 - X2)T exp li d;z°e-*^° V(r)^(x° - z'^)e'''^ 
xe->l{e{x[ - xf) + ^yi^i - xi)7ioe^<^' 



xTexp l-^J^ dz''e-''°'''V{r')e{z^ - x'^'^^'j. (18) 
Using for the delta-function the closure relation: 

<53(xi-x;)=5:^„(x)^„(x') (19) 

n 

[from the expression of G20, Eq. (^), we have X2 = xy, where the tp^s form a complete 
set of spinor wave functions, and comparing the resulting expression of G with its cluster 
decomposition formula (for and x^ greater than xf and x'^) |1T7| : 

G(xi, X2, x[,x'2) = $„(xi, X2)$„(x;, x'2) = 0n(x)e-'^-(^-^') 0„(x'), (20) 

n n 

where the $'s are generalized Bethe-Salpeter wave functions, one deduces the energy 
spectrum of the theory. [We have introduced the total momentum variable P = pi + P2 
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and the total coordinates X = {xi + X2)/2 and X'.] By identifying in the expression of G 
the argument of the exponential functions containing the total coordinates and X'^ 
with the energy eigenvalues, one deduces that the energy spectrum is that of the positive 
energy solutions of the Dirac equation with the static potential V (added by 1712): 

(^7i.9i - mi - 7ioV(r)) VK, xi - xa) = 0. (21) 

Repeating the above operations with values of X2 and xf greater than X2 and x^, one 
recovers the negative energy solutions of Eq. (0). Therefore, the energy spectrum of 
the two-particle Green's function (^ in the limit 777,2 — ^ C)0 is given, up to the additive 
factor 7772, by the spectrum of eigenvalues of Eq. (^). This spectrum does not contain 
any abnormal solutions corresponding to the relative time excitations. 

The previous procedure makes also possible the identification of the Bethe-Salpeter 
wave function in terms of the Dirac wave function ip. Designating by piQ the energy 
eigenvalue of the Dirac wave function, one has for the bound state wave functions (for 
which Xi > x'l): 

$(xi,X2) = e-^("2^2+Pi«^?)Texp i^i £ dz^e-''°^V{r)e{x'^ - z°)e*"°'^|^(x). (22) 
For x^ < 0, the integral vanishes and one obtains: 

$ _ ^-i(m2X°2+P^0X°)^(^^^ ^ g-i(mi+m2+E)X''g-i(mi-m2)xO/2g-iEx''/2^^^^^ ^0 ^ ^23) 

where we have introduced the binding energy E: 

Pio = mi + E, E <Q. (24) 
For x^ > 0, the ^-function in the integral can be replaced by 1, and using the relation p8 



^{A+B)t^-At 



Texp^*(it'e^*'5e-^*'|, (25) 



one finds: 



$ = g-*(™2+Pio)a:^g-ia;"(/i-y)^j^^'j 

= ^-i(mi+m2+E)X" ^-i(mi-m2)x° /2 ^iEx° /2 ^-ix°{ho-mi) ^^^^^ ^ ^0 ^ j^2g) 

In the bound state domain, the operator {Hq — mi) has a positive spectrum since it is 
essentially equivalent to the kinetic energy operator. One can then check that the above 
Bethe-Salpeter wave function has the correct spectral properties 0, H. For x° > 0, 
it has singularities for positive values of the relative energy variable, and for x° < 0, 
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it has singularities for negative values of the relative energy variable. Furthermore, by 
decomposing $ as 

$ = ^(a;0)$+ + ^(-x°)<l>_, (27) 

one verifies that the conjugate Bethe-Salpeter wave function $, identified from the com- 
parison of Eqs. (0) and (20) with the use of Eq. ([T9|) , satisfies the correct definition 
(with the antichronological product): 



$ = ^(-xO)$t7io + ^(x")$L7io. 



(28) 



The Bethe-Salpeter wave function could also have been constructed from the resolu- 
tion of the equations it satisfies. By taking in Eq. (^ the limit X2 — — 00 and using the 
cluster decomposition (pO]), one finds the equations: 



i^i.di - mi - 7ioV(r)^(x° - x?))<l> 



0, 



(29) 



the solution of which is precisely the function defined in Eqs. ( |23|) and (pq) , when the 
continuity condition at x° = is used. 
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4 Covariant gauges 



For completeness, we briefly sketch the derivation in the case of hnear covariant gauges. 
Here, a comphcation arises because of the presence of gauge photons (scalar and longitu- 
dinal) which create in the spectrum of states additional cuts starting from the bound state 
poles |TP[ . These gauge dependent cuts are usually removed by appropriately associating, 
through Ward-Takahashi identities, with the exchanged photons the photons contribut- 
ing to the radiative corrections and taking in the cluster decomposition formula ( pOt ) the 
limits CX3 and X'^ —>■ — oo, while keeping and finite. With this procedure one 

selects the bound state spectrum |T^, which is the object of interest. An illustration of 
this cancellation can be found in Ref. |jl2| for the static case. We shall not consider here 
the radiative corrections, which are more complicated than in the static limit, but shall 
indicate when necessary, their precise role. 

In linear covariant gauges, characterized by a gauge parameter ^, the photon propa- 
gator in momentum space is: 

C..(*) = -(ft..-f^)^. (1) 

Defining the functions 

x'{x) = -f--^ In ( r: ' ^„ : r i , (2) 




[x' being the derivative, in the Feynman gauge, with respect to x^, of the function x 
defined in Ref. [|l3;] one finds for the effective potential [Eq. (0)] the expression: 



A\^\x,X2,X2) = ig ^,o(x'ix - X2) - ie^^d^(^A{x - X2) - A{x - x'^)^. (4) 



The fact that the ^-dependent part od A^^^ is a derivative term, makes it possible its 
immediate elimination from Eqs. One finds for G^^^ the following ^ dependence: 

G^^) = exp i^ie^^(^A{xi - X2) - A(xi - x'^) - A{x[ - X2) + A{x[ - x'^)^^G^°\ (5) 

which is precisely the general transformation law of the two-particle Green's function 
under covariant gauge transformations of the photon propagator (radiative corrections 
being neglected here) [^, It has been shown that these ^-dependent gauge 

factors do not change the bound state pole positions of Green's functions [|T5|, [T6|. It is 
therefore sufficient to work in the Feynman gauge (^ = 0), which is also similar to the 
scalar interaction case. (We shall henceforth omit the index from the Green's function.) 
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The function x' has the following asymptotic behavior: 

1 „ , 1 



X'{x) ^ e{x')-V{\^\) + 0{-). (6) 



2 

It is the 0(1/ x^) term which is gauge dependent (it is proportional to the derivative of 
the function A of Eq. (0) and vanishes in the gauge ^ = —2) and must be cancelled by 
a corresponding piece coming the radiative corrections [1^. We shall assume that this 
cancellation has occurred and thus the next-to-leading term in the asymptotic part of x' 
is O(l/a;02). 

Replacing A^^^ by its expression (|^) in Eqs. (||)-(H) and isolating again the asymptotic 
part of the potential, as in Eqs. (|1^-(|15]), one obtains: 



(^i'ji.di - mi - 7ioV'(r) - 710 (^x'{x2 - Xi) - ^V{r)^ + (^xi^i - x'^) - ^V{r) 



G 



= G2o{x2 - X2)i6 {xi - x[), (7) 

and a similar conjugate equation. The remaining part of the resolution of these equations 
parallels that utilized in the Coulomb gauge. One arrives at an expression analogous to 
that of Eq. (0), where the potential part V{1 — 99) is replaced by the bracket term 
of Eq. (with x\ replaced by the integration variable y\). One then takes the limits 
x'j' ^ 00 and X2-, xf —00, by keeping x^ and finite. Because the difference 
x'{x) — V/2 behaves as 0(l/x°^), in the part containing x'(x2 — yi) — V/2 it is the values 
of ?/i close to X2 that contribute, while in the part containing x'ivi — 3:2) — V^/2 it is the 
values of close to X2 that are relevant. Thus these two parts receive contributions from 
disjoint intervals and the T-product factorizes again [Eq. (plTf )], with the first integral 
lying from —00 to Xi and the second integral from xf to 00. One ends up, for the bound 
state spectrum, with the same conclusion as in the Coulomb gauge. 
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